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INTRODUCTION 



Many of you are vigorously searching for the quark-gluon plasma — a predicted 
new phase of nuclear matter where quarks roam almost freely throughout the medium 
instead of being confined to individual nucleons. 1 ' 2 Such a plasma is believed to have 
existed in the first 10 /is of the universe during the big bang and could be produced in 
the laboratory during the little bang of a relativistic heavy-ion collision. 

When nuclei collide head-on at relativistic speeds, the nuclear matter is initially 
compressed and excited from normal nuclear density and zero temperature to some 
maximum values — during which pions, kaons, and other particles are produced — and 
then expands, with a decrease in density and temperature. The early stages of the 
process are often treated in terms of nuclear fluid dynamics, but at some late stage the 
expanding matter freezes out into a collection of noninteracting hadrons. 

To sample the density, temperature, collective velocity, size, and other properties 
of the system during this freeze-out, some of you are measuring invariant one-particle 
multiplicity distributions and two-particle correlations for the pions, kaons, and other 
particles that are produced. Your hope is that a sharp discontinuity in the value 
of one or more of the extracted freeze-out properties as a function of bombarding 
energy and/or size of the colliding nuclei could signal the formation of a quark-gluon 
plasma or other new physics. For the extraction of these freeze-out properties from your 
experimental measurements, a nine-parameter expanding source model was presented 
at the 12th Winter Worshop on Nuclear Dynamics. 3 



NINE-PARAMETER EXPANDING SOURCE MODEL 



This source model describes invariant one-particle multiplicity distributions and 
two-particle correlations in nearly central relativistic heavy-ion collisions in terms of 
nine parameters, which are necessary and sufficient to characterize the gross properties 
of the source during its freeze-out from a nuclear fluid into a collection of noninter- 
acting, free-streaming hadrons. 3-5 The values of these nine parameters, along with 
their uncertainties at 99% confidence limits, are determined by minimizing x 2 for the 
types of data considered. Several additional physically relevant quantities, along with 
their uncertainties at 99% confidence limits, can then be directly calculated. The nine 
independent source freeze-out properties that we consider here are the central baryon 
density n, nuclear temperature T, transverse collective velocity v t , longitudinal collec- 
tive velocity v e , source velocity v s , transverse radius i? t , longitudinal proper time Tf, 
width in proper time At, and pion incoherence fraction A„-. 

For a particular type of particle, the invariant one-particle multiplicity distribution 
and two-particle correlation function are calculated in terms of a Wigner distribution 
function S(x,p), which is the phase-space density on the freeze-out hypersurface, giving 
the probability of producing a particle at spacetime point x with four-momentum p. 
It includes both a direct term 6 and a term corresponding to 10 resonance decays. 7 
We consider nearly central collisions, assume axial symmetry, and work in cylindrical 
coordinates in the source frame, with longitudinal distance denoted by z, transverse 
distance denoted by p, and time denoted by t. Throughout the paper we use units in 
which h — c — k — 1, where H is Planck's constant divided by 2n, c is the speed of 
light, and k is the Boltzmann constant. However, for clarity, we reinsert c in the units 
of quantities whose values are given in the text or table. 

Integration of the direct term over spacetime leads to the Cooper-Frye formula for 
the direct contribution to the invariant one-particle multiplicity distribution, 8 namely 

p d . ( ) = E ^dir = 1 ^dir = 2J+1 I ^ p» 

dp 3 2ixm t dydm t (2ir) 3 7s ^ exp{[p • v(x) — p(x)]/T(x)} =p 1 

(1) 

where E denotes the particle's energy, m t = \Jw? + p t 2 its transverse mass, and y 
its rapidity. The quantity m is the particle's rest mass, and p t = \Jp x 2 + P y 2 is its 
transverse momentum. The minus sign applies to bosons and the plus sign to fermions. 
The quantity J is the spin of the particle, v (x) is the collective four-velocity, T(x) is 
the nuclear temperature, and p(x) is the chemical potential for this type of particle. 
We assume that the source is boost invariant within the limited region between its two 
ends, 9 ' 10 and that it starts expanding from an infinitesimally thin disk at time t = 0. 
The transverse velocity at any point on the freeze-out hypersurface, whose integration 
limits are denoted by S, is assumed to be linear in the transverse coordinate p. 

For a particular type of particle, the two-particle correlation function is given 
by 4 ' 5,11 

C(K a\ p 2(pi,p 2 ) _ 1+ A \I d 4 xS(x,K) expjiq- x)\ 2 
{ > q) P( Pl )P( P2 ) [fd*xS(x, Pl )][fd*xS(x,p 2 )}' {) 

where K = \{j>\ + P2) is one-half the pair four-momentum and q = pi — p 2 is the pair 
four-momentum difference. In this equation the plus sign applies to bosons and the 
minus sign to fermions, and the quantity A specifies the fraction of particles of this type 
that are produced incoherently. 
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APPLICATION TO Pb + Pb COLLISIONS AT p lah /A = 158 GeV/c 

We have used the nine-parameter expanding source model described in the pre- 
vious section to analyze normalized but still preliminary data from Experiment NA44 
performed at CERN's Super Proton Synchrotron. 12 ' 13 These data consist of invariant 
7r + , K + , and K~ one-particle multiplicity distributions and 7r + and n~ two-particle cor- 
relations for the central 20% collisions in the reaction Pb + Pb at piab/A = 158 GeV/c. 

For symmetric collisions the source velocity v s can be calculated in terms of the 
beam momentum per nucleon and nucleon mass, which eliminates the need to vary this 
parameter. The remaining eight adjustable parameters are determined by minimizing 
X 2 with a total of 2137 data points for the five types of data considered, so the number of 
degrees of freedom v is 2129. The error for each point is calculated as the square root of 
the sum of the squares of its statistical error and its systematic error, with a systematic 
error of 15% for the 7r + , K + , and K~ one-particle multiplicity distributions and zero 
for the 7r + and 7i~ two-particle correlations. The resulting value of x 2 is 2165.5, which 
corresponds to an acceptable value of x 2 l u — 1-017. There is a 28.6% probability that 
X 2 '/v would be at least this large for a perfect model. The individual values of x 2 l v are 
0.904, 0.925, and 0.721 for the 7r + , K + , and K~ one-particle multiplicity distributions 
and are 1.130 and 1.012 for the 7r + and 7r~ two-particle correlations, respectively. 

The values of the independent freeze-out properties determined this way, along 
with their uncertainties at 99% confidence limits on all quantities considered jointly, 
are given in the third column of Table 1. For comparison, we show in the second column 
of Table 1 the corresponding results 4 for the central 7% collisions in the reaction Si + Au 
at = 14.6 GeV/c studied 14 ' 15 in Experiment E-802 at Brookhaven's Alternating 

Gradient Synchrotron. Compared to these earlier results for Si + Au collisions, the 
present results indicate that in Pb + Pb collisions the freeze-out density is somewhat 
lower, the freeze-out temperature is slightly higher, the source at freeze-out is somewhat 
larger, and the longitudinal collective velocity is very poorly determined (because of 
the limited experimental coverage in rapidity). The quantity n appearing in Table 1 
denotes normal nuclear density, whose value 16 is calculated from the nuclear radius 
constant r by means of n = 3/(47rr 3 ) = 3/[47r(1.16 fm) 3 ] = 0.153 fm~ 3 . 

Table 1. Comparison of nine independent source freeze-out properties for 
the central 20% collisions of Pb + Pb at pub/A = 158 GeV/c with those 
for the central 7% collisions of Si + Au at piab/A = 14.6 GeV/c. 



Value and uncertainty at 99% confidence 
Property Si + Au Pb + Pb 



Central baryon density n/no 


U.140 _o.o4 5 


062 + 019 

U.UOZ _ .015 


Nuclear temperature T (MeV ) 


92.9 ± 4.4 


95.8 ± 3.5 


Transverse collective velocity v t (c) 


0.683 ± 0.048 


0.664 ± 0.035 


Longitudinal collective velocity vg (c) 


900 + 023 
u.yuu _ .029 


0.9985 1°;°4 15 


Source velocity v s (c) 


875 + - 015 

U.O/O _ 016 


0.9941219 (fixed) 


Transverse radius R t (fm) 


8.0 ± f.6 


11.4 ± f.5 


Longitudinal proper time Tf (fm/c) 


8.2 ± 2.2 


12.2 ± 2.1 


Width in proper time At (fm/c) 


5 q +4-4 


7 1 + 4 - 5 

1.1 -2.4 


Pion incoherence fraction X v 


0.65 ± O.ff 


0.690 ± 0.074 
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Figure 1. Comparison between model predictions and experimental data for the invariant 
one-particle multiplicity distribution Ed 3 N/dp 3 = l/(27rm t ) d 2 N / 'dy drat for n + , K + . and K~ as a 
function of rat — fn for y = 2.675. The error bars shown in this figure represent statistical errors only. 



Figure [T] shows an example of our model predictions, given by the curves, compared 
with experimental data for the invariant one-particle multiplicity distribution for tt + , 
K + , and K~ . Solid lines and solid symbols are used for positive particles, and a dashed 
line and open symbols are used for negative particles. The three curves are calculated 
for a value of rapidity y = 2.675, which is the central value of the kaon rapidity coverage 
in this experiment. The experimental data for kaons correspond to this same rapidity, 
with a rapidity bin width Ay = 0.05. For pions, the high-m t experimental data (those 
to the right of the noticeable break in the distribution of points) correspond to this 
same rapidity and bin width, whereas the low-mt experimental data correspond to the 
nearby rapidity y = 2.65 and rapidity bin width Ay = 0.10. The data shown in this 
figure represent only a small fraction of that used in our analysis. 

Since the two-particle correlation function depends on five variables, it is some- 
what more difficult to graphically compare our model predictions with experimental 
data. The comparisons that we have made thus far involve fixing the pair rapidity 
and transverse momentum to specific values and then plotting the two-particle corre- 
lation function versus each component of the pair three-momentum difference, namely 
longitudinal, Qside, and g out , for fixed values of the other two components. These com- 
parisons demonstrate that our expanding model satisfactorily reproduces experimen- 
tal two-particle correlations. However, as indicated earlier, the agreement between 
model predictions and experimental data is somewhat better for negative pions, where 
X 2 /v = 1-012, than for positive pions, where \ 2 l v = 1-130. 
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RECONCILIATION WITH PREVIOUS ANALYSES 



Analyses with our expanding source model for both the reaction Si + Au at 
Pi&b/A = 14.6 GeV/c and the reaction Pb + Pb at piab/A = 158 GeV/c indicate that 
the freeze-out temperature is less than 100 MeV and that both the longitudinal and 
transverse collective velocities — which are anti-correlated with the temperature — are 
substantial. Similar conclusions concerning a low freeze-out temperature have also been 
reached by Csorgo and Lorstad. 17,18 However, other analyses 12,13,19-21 have yielded a 
much higher freeze-out temperature of approximately 140 MeV. In order to reconcile 
this serious discrepancy, we now examine the features in these analyses that led them 
to the conclusion of a much higher freeze-out temperature. These analyses fall into two 
major classes, which we consider in turn. 

Neglect of Relativity in Extrapolation of Slope Parameters 

One type of analysis 12,13 was based upon the extrapolation to zero particle mass of 
extracted slope parameters characterizing the dependence of unnormalized transverse 
one-particle multiplicity distributions upon transverse mass. For a given reaction and 
type of particle, this transverse one-particle multiplicity distribution was represented 
by the expression?] 

1 dN . ( m t 

where A is an arbitrary normalization constant and T e g is the extracted slope parameter. 
Values of T e g were extracted in this way for six types of particles originating from three 
separate reactions, namely tt + , 7r~, K + , K~, p, and p originating from the reaction 
p + p at center-of-mass energy y/s = 23 GeV, from the central 10% collisions in the 
reaction S + S at p^b/A = 200 GeV/c, and from the central 6.4% collisions in the 
reaction Pb + Pb at p\ a b/A = 158 GeV/c. 

As we will see below, the values of these extracted slope parameters contain valu- 
able information, but they were unfortunately analyzed in Refs. 12 and 13 in terms of 
a heuristic equation that neglects relativity, namely 

T eff = T + mv 2 , (4) 

where T is the nuclear temperature (whose value we are trying to determine) and 
v is the average transverse collective velocity of the expanding matter from which the 
particle originated. On the basis of Eq. (Q), the extrapolation in Ref. 13 of the extracted 
slope parameters to zero particle mass yielded the result T w 140 ± 15 MeV. 

In the limit in which the particle velocity is large compared to the average col- 
lective velocity and with the aid of other simplifying assumptions and approximations, 
the correct relationship between slope parameter, nuclear temperature, particle mass, 
and average collective velocity can be easily derived from the relativistically correct 
Eq. ([!]). With the neglect of contributions from resonance decays, the neglect of the 
=Fl appearing in the denominator of Eq. ([!]), the assumption of a constant freeze-out 
temperature, and the assumption that freeze-out occurs at a constant time t in the 

*To facilitate comparisons with our own expressions, we have transformed the notation used in 
Refs. 12 and 13 to that used here. 
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source frame, Eq. (|l|) leads to 

7 (x)[£-p.v(x)] 



m t dydm t 



/ d 3 x exp 


p ■ u(x) 


= A'E J eft; exp j 


IV 


T 





T 



(5) 

where A' is a different arbitrary normalization constant from the one appearing in 
Eq. (^), the subscript V on the integral denotes the spatial integration limits for the 
source, and the position-dependent Lorentz factor 7(x) = 1/jT— v(x)-v(x). 

By introducing an average collective velocity v in the integrations in Eq. (||), taking 
the limit in which the particle velocity is large compared to the collective velocity, spe- 
cializing to the transverse direction, and neglecting the pre-exponential E dependence, 
we are led to 

7(m t - p t v) 



1 dN 

: — = A exp 

m t dm t 



T 



( m t - v^rrit 2 - m 2 \ 

Aexp { IVT - P j ' (6) 



To obtain the relationship between the slope parameter, nuclear temperature, particle 
mass, and average collective velocity, we equate the derivatives with respect to m t of 
Eqs. (H) and which leads to 

T=(l-^)^=(l-,^)^=. (7) 



Pt J vi - v 2 V V Vt 2 I \f\ 



v 



An analogous relationship has also been obtained by Siemens and Rasmussen 22 for the 
case of a blast wave produced by the explosion of a spherically symmetric fireball. 
In the limit of zero particle mass, Eq. (0) reduces to 



'-eS\ 



'1 - V 

1 + v 



which agrees with the result obtained by Schnedermann, Sollfrank, and Heinz 23,24 for 
the case of cylindrical symmetry. With a typical value of 0.4 c for the average collective 
velocity v and the limiting value of T e fj w 140 ± 15 MeV obtained in Ref. 13 by 
extrapolating slope parameters to zero particle mass, Eq. @ yields T w 92 ± 10 MeV 
for the nuclear temperature. 

Several Approximations Made in a Thermal Model 

Another type of analysis 12,13,19-21 utilized the thermal model of Schnedermann, 
Sollfrank, and Heinz 23,24 to extract the nuclear temperature and transverse surface 
collective velocity from unnormalized experimental transverse one-particle multiplicity 
distributions. An accumulation of effects from several approximations led to a some- 
what higher temperature than we have found with our expanding source model. These 
approximations include the neglect of contributions from resonance decays, the neglect 
of the =Fl appearing in the denominator of Eq. ([!]), the neglect of the coupling of the 
transverse motion to the longitudinal motion, and — most importantly — the neglect of 
information contained in the absolute normalization of the multiplicity distributions. 
The accumulation of effects from these approximations was responsible for the conclu- 
sion on page 2083 of Ref. 13 that "Within a temperature range 100 <T< 150 MeV, the 
fits are equally good." It is seen that the use of unnormalized experimental transverse 
one-particle multiplicity distributions in such a thermal model provides only a rough 
indication of the nuclear temperature at freeze-out. 
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ADDITIONAL STUDIES WITH EXPANDING SOURCE MODEL 



To test the robustness of our expanding source model, we used it to analyze one- 
particle and correlation data generated theoretically from a nuclear fluid dynamical 
calculation performed with the computer program HYLANDER 25 that corresponded to a 
high freeze-out temperature and an extremely low transverse freeze-out velocity. The 
results of this study demonstrated that the model is capable of reproducing the under- 
lying freeze-out properties even when their values were chosen to lie in unanticipated 
regions. 

To determine whether or not ultrarelativistic proton-proton collisions can be de- 
scribed in terms of nuclear fluid dynamics, we used our expanding source model to 
analyze invariant 7r + , 7r~, K + , K~ , p, and p one-particle multiplicity distributions 26 for 
the reaction p + p at center-of-mass energy y/s = 45 GeV. We included the systematic 
and normalization errors discussed in Ref. 26 in addition to statistical errors. Because 
the pion incoherence fraction does not enter in the expression for one-particle multi- 
plicity distributions and because for symmetric collisions the source velocity v s can be 
calculated, there are only seven adjustable parameters in this case. These parameters 
are determined by minimizing x 2 with a total of 459 data points for the six types of 
data considered, so the number of degrees of freedom v is 452. The resulting value of 
X 2 is 1132.0, which corresponds to a completely unacceptable value of x 2 l v — 2.504. 
The probability that a perfect model would have resulted in a value of x 2 & t least as 
large as that found here is the incredibly small value 4.9 x 10~ 60 . 



SUMMARY AND CONCLUSIONS 

We have used a nine-parameter expanding source model that includes special rel- 
ativity, quantum statistics, resonance decays, and freeze-out on a realistic hypersurface 
in spacetime to analyze in detail invariant 7r + , K + , and K~ one-particle multiplicity 
distributions and tt + and 7r~ two-particle correlations in nearly central collisions of 
Pb + Pb at pi a b/^4 = 158 GeV/c. These studies confirm an earlier conclusion; 3-5 for 
nearly central collisions of Si + Au at pi a b/A = 14.6 GeV/c the freeze-out temperature is 
less than 100 MeV and both the longitudinal and transverse collective velocities — which 
are anti-correlated with the temperature — are substantial. 

We also reconciled our current results with those of previous analyses that yielded 
a much higher freeze-out temperature of approximately 140 MeV for both Pb + Pb 
collisions at p\ a h/A = 158 GeV/c and other reactions. One type of analysis was based 
upon the use of a heuristic equation that neglects relativity to extrapolate slope pa- 
rameters to zero particle mass. Another type of analysis utilized a thermal model in 
which there was an accumulation of effects from several approximations. 

The future should witness the arrival of much new data on invariant one-particle 
multiplicity distributions and two-particle correlations as functions of bombarding en- 
ergy and/or size of the colliding nuclei. The proper analysis of these data in terms 
of a realistic model could yield accurate values for the density, temperature, collective 
velocity, size, and other properties of the expanding matter as it freezes out into a 
collection of noninteracting hadrons. A sharp discontinuity in the value of one or more 
of these properties could conceivably be the long-awaited signal for the formation of a 
quark-gluon plasma or other new physics. 
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